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Abstract 

Data  for  a  scene  of  interest  may  be  collected  over  multiple  polarization  channels. 
In  the  case  of  polarimetric  synthetic  aperture  radar  images,  regularization  techniques 
are  typically  applied  independently  to  each  polarimetric  channel.  However,  independent 
processing  does  not  account  for  cross-channel  coupling  and  may  corrupt  the  polarimetric 
information  in  the  signals.  Recent  consideration  of  joint  enhancement  techniques 
has  shown  promising  results  for  multi-channel  datasets  with  similar  regions  of  signal 
magnitude  and/or  phase.  However,  in  the  case  of  polarimetric  SAR  data,  scattering  may  be 
present  in  some  channels  and  not  in  others.  This  thesis  mathematically  formulates  multi¬ 
channel  sparse  imaging  for  polarimetric  radar  data  using  a  joint  enhancement  algorithm  to 
enforce  sparsity  and  polarimetric  coupling  constraints. 

Two  candidate  functional  relationships  are  derived  to  describe  polarimetric  coupling 
among  received  signal  channels:  one  convex  function  and  one  non-convex  function.  These 
functions  are  reformed  as  optimization  constraints.  Then,  an  optimization  problem  is 
constructed  to  maintain  signal  fidelity,  enforce  sparsity,  and  preserve  interchannel  coupling. 
An  iterative  dual  gradient  descent  algorithm  is  used  to  alternatively  calculate  updated  scene 
estimates  for  each  channel  and  the  maximizing  Lagrange  multipliers  for  each  coupling 
constraint.  Results  are  found  for  several  polarimetric  SAR  datasets.  Jointly  enhanced 
images  are  compared  with  corresponding  images  found  through  independent  enhancement, 
taking  into  consideration  signal  fidelity,  sparsity,  polarimetric  preservation,  and  scattering 
classification. 

Overall,  the  jointly  enhanced  image  channels  display  significantly  better  polarimetric 
preservation  compared  to  the  corresponding  independently  restored  image  channels.  More 
research  is  needed  to  understand  how  improved  polarimetric  preservation  can  be  used  to 
improve  target  classification. 
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ENHANCED  POLARIMETRIC  RADAR  IMAGING  USING  CROSS-CHANNEL 


COUPLING  CONSTRAINTS 

I.  Introduction 

IN  the  late  1800’s,  physicist  Heinrich  Hertz  demonstrated  several  important  properties  of 
electromagnetism  that  sparked  widespread  research  interest  in  radio  waves.  Among 
these  experiments,  Hertz  showed  that  radio  waves  can  be  transmitted,  re-radiated  by 
metallic  objects,  and  then  detected  [1].  These  important  discoveries  led  to  the  development 
of  Radio  Detection  and  Ranging  (RADAR).  Early  radar  systems  were  simple  and  primarily 
used  to  determine  a  target’s  range  and  velocity  in  a  given  direction  [2].  Modem  radar 
systems  also  have  the  capability  to  track  targets,  determine  a  target’s  size  and  shape,  form 
images  of  a  scene  of  interest,  and  eliminate  unwanted  interference.  Principles  of  radar 
have  led  to  advancements  in  many  fields  including  meteorology,  law  enforcement,  military 
operations,  and  medicine  [2]. 

The  importance  of  wave  polarization  in  radar  was  first  studied  extensively  by 
George  Sinclair  in  the  1940’s  [1].  Sinclair  determined  that  radar  targets  transform  the 
transmitted  wave  polarization  based  on  target  attributes  such  as  material  and  shape.  By 
transmitting  and  receiving  both  a  Horizontal  (H)  and  Vertical  (V)  polarization,  Sinclair 
characterized  basic  properties  of  targets  of  interest  [1].  Sinclair’s  foundational  work  led 
to  significant  advances  in  polarimetry,  the  study  of  target  aspects  that  can  be  determined 
by  changing  the  transmitted  and  received  wave  polarization  state  [1].  State-of-the-art 
radar  systems  use  multiple  transmit/receive  polarization  channels  for  automated  target 
classification  and  interference  reduction  [3].  This  thesis  presents  a  method  to  preserve 
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the  relationship  between  multiple  radar  polarization  channels  during  sparsity-constrained 
image  enhancement  in  order  to  improve  scene  estimates  and  target  classification. 

1.1  Motivation 

Synthetic  aperture  radar  (SAR)  images  contain  impulse  response  effects  that  appear  to 
smear  or  blur  the  image  [4].  Sparse  imaging  is  a  regularization  method  which  has  proved 
desirable  for  many  radar  applications  since  it  provides  a  way  to  remove  smearing  resulting 
from  limited  datasets  [3].  In  essence,  regularization  methods  estimate  scene  reflectivity  by 
reducing  the  impulse  response  effects.  Sparsity  constraints  enforced  alongside  data  fidelity 
constraints  lead  to  signal  energy  placement  in  only  high  amplitude  pixels  [3].  Hence, 
enforcing  sparsity  in  an  image  emphasizes  high  magnitude  scatterers  while  reducing  the 
effects  of  noise  and  the  impulse  response. 

Often  times,  data  for  a  scene  of  interest  is  collected  over  multiple  channels,  such 
as  frequency,  elevation,  or  polarization  [5].  In  the  case  of  polarimetric  SAR  images,  the 
regularization  techniques  are  typically  applied  independently  to  each  polarimetric  channel. 
However,  independent  processing  may  corrupt  the  polarimetric  information  in  the  signals, 
as  it  does  not  account  for  cross-channel  coupling.  Recent  research  into  joint  enhancement 
techniques  has  shown  promising  results  for  multi-channel  datasets  with  similar  support 
regions,  i.e.,  regions  of  similar  signal  magnitude  and/or  phase  [5], [6].  However,  in  the  case 
of  polarimetric  SAR  data,  scattering  may  be  present  in  some  channels  and  not  in  others. 
This  thesis  mathematically  formulates  multi-channel  sparse  imaging  for  polarimetric  radar 
data  using  a  joint  enhancement  algorithm  to  enforce  sparsity  and  polarimetric  coupling 
constraints. 

1.2  Layout 

Chapter  2  presents  relevant  background  information  on  polarimetry,  SAR,  and  current 
image  enhancement  techniques.  The  Sinclair  scattering  matrix  is  presented  along  with 
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related  polarimetric  decompositions,  which  are  later  used  to  formulate  coupling  constraints. 
Moreover,  current  techniques  for  independent  and  joint  image  enhancement  are  presented, 
as  well  as  the  problems  that  arise  with  each  when  considering  polarimetric  data. 

Chapter  3  derives  two  possible  functional  relationships  to  describe  polarimetric 
channel  coupling  among  received  signal  channels.  Then,  an  optimization  problem  is 
constructed  to  maintain  signal  fidelity,  enforce  sparsity,  and  preserve  interchannel  coupling. 
The  convexity  of  the  problem  is  analyzed.  Finally,  an  algorithmic  solution  of  the 
optimization  problem  is  developed  and  a  measure  of  scattering  classification  is  defined. 

In  Chapter  4,  the  joint  optimization  problem  is  solved  using  the  algorithm  described 
in  Ch.  3  for  several  polarimetric  SAR  datasets.  These  datasets  include  simulated  radar  data 
of  various  canonical  scatterers.  Jointly  enhanced  images  are  compared  with  corresponding 
images  found  through  independent  enhancement,  taking  into  consideration  signal  fidelity, 
sparsity,  polarimetric  preservation,  and  scattering  classification. 

Chapter  5  concludes  the  findings  of  our  research  and  provides  ideas  for  future  research 
into  this  topic. 
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II.  Background 


Tnis  chapter  introduces  the  basics  of  polarimetry,  synthetic  aperture  radar  (SAR),  and 
sparse  imaging.  Then,  the  multichannel  imaging  problem  is  introduced,  including 
current  methods  for  independent  and  joint  image  enhancement. 

2.1  Notation 

Scalars  are  denoted  by  italicized  letters.  Vectors  and  matrices  are  bold  and  single  and 
double  underlined,  respectively,  such  as  /  and  M.  The  /th  entry  in  a  vector  is  given  by  /(/). 
Similarly,  the  i yth  entry  in  a  matrix  is  denoted  M(i,  j).  A  set  of  n  vectors  or  matrices  is 
denoted  by  a  subscript  such  as  /  ,  /  ,...,/  .  The  term  diag  {•}  represents  a  diagonal  matrix 
whose  entry  is  defined  by  the  expression  in  the  brackets.  The  operator  (-j^  represents 
the  Hermitian  operator  for  vectors  and  matrices,  and  represents  the  complex  conjugate  for 
scalars. 


2.2  Polarimetry 

In  a  monostatic  radar  system,  the  transmitting  and  receiving  antennas  are  colocated 
[2].  An  electromagnetic  wave  pulse  is  transmitted,  interacts  with  a  target  of  interest,  and 
is  reradiated  non-uniformly  in  every  direction.  The  backscattered  pulse  is  collected  at  the 
receiver.  The  ideal  signal  power  received  at  the  antenna,  Pr,  is  given  by  the  radar  range 
equation 


PtGWcr 


(2.1) 


where  Pt  is  the  transmitted  signal  power,  G  is  the  antenna  gain,  w  is  the  wavelength,  cr  is 
the  target  radar  cross  section  (RCS)  and  R  is  the  one-way  range  from  the  radar  to  the  target 
[2].  The  RCS  is  determined  by  target  characteristics  such  as  size,  shape,  or  material,  and 
radar  parameters  such  as  wavelength  and  incidence  angle  [7].  The  RCS  is  a  measure  of 
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detectability  and  hence  provides  an  effective  area  of  the  target  [7,  8].  Definitively, 

cr  =  lim  «  4;7r|5  p  in  the  far  field,  (2.2) 

R^oo 

where  represents  the  intensity  of  the  transmitted  and  received  electromagnetic  waves, 
respectively,  and  |5|  is  defined  as  the  complex  scattering  amplitude  of  the  wave  [7,  8]. 
2.2.1  Scattering  Matrix. 

The  RCS  cr  also  depends  on  the  polarization  of  the  incident  and  scattered  waves. 
Specifically, 

(Tpq  «  4;r|S^<^P,  (2.3) 

where  p  is  the  polarization  of  the  transmitted/incident  field  and  q  is  the  polarization  of  the 
scattered/received  field  [1][7].  The  polarization  states  p  and  q  are  orthogonal,  such  as  a 
linear  basis  with  horizontal  and  vertical  components  or  circular  basis  with  left  and  right 
components.  For  simplicity,  this  work  assumes  p,  q  are  either  horizontal  (H)  or  vertical 
(V).  The  transmitted  and  received  fields  have  components  in  the  horizontal  direction  h  and 
vertical  direction  v,  i.e., 

=  Elv  +  E]jh,  (2.4) 

E^  =  Ell  +  Elh,  (2.5) 

where  Eh,  Ey  denote  the  horizontal  and  vertical  components,  respectively,  of  the  incident 
and  scattered  fields  [8].  A  fully  polarimetric  radar  can  transmit  and  receive  both  horizontal 
and  vertical  polarizations.  By  alternating  the  transmit  and  receive  polarization  states,  the 
radar  collects  all  four  complex  scattering  amplitudes.  The  entire  scattering  response  is 
stored  in  the  Sinclair  or  Scattering  Matrix  S,  which  relates  the  transmitted  and  received 
fields  as  #  =  [8]  [9].  Explicitly, 


^ikoR 

s  HH 

s  HV 

E^ 

V  y\ 

R 

S  VH 

S  vv 

V  v\ 
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where  Uq  is  the  wavenumber  of  the  transmitted  wave,  and  the  term  accounts  for 
magnitude  and  phase  deviations  due  to  propagation  over  the  distance  R  [9]. 

The  four  scattering  amplitudes  in  5  are  complex-valued  and  hence  contain  both 
a  magnitude  and  relative  phase.  Each  scattering  amplitude  is  a  function  of  the  radar 
frequency  and  the  target  aspect  angle  [8].  The  diagonal  elements  are  known  as  co¬ 
polarization  terms  since  they  relate  the  scattering  amplitude  for  the  same  transmit  and 
receive  polarization.  Similarly,  the  olf-diagonal  elements  are  known  as  cross-polarization 
terms  since  they  relate  perpendicular  polarizations.  In  the  case  of  monostatic  radar,  the 
cross-polarization  terms  are  equal  due  to  electromagnetic  reciprocity  [7-9].  Note  that  a 
backscatter  coordinate  system  is  assumed,  which  is  radar-centric.  Within  this  model,  the 
horizontal  and  vertical  transmit  channels  are  initially  olf-set  in  phase  by  180°  or  n  radians 
[8,  9]. 


The  material  and  geometry  of  a  particular  target  determine  the  characteristics  of  the 
reradiated  wave.  When  an  incident  pulse  reflects  once  off  a  target  surface,  the  signal  is 
said  to  single  bounce.  Similarly,  when  an  incident  pulse  is  scattered  by  two  surfaces  of 
the  target  before  returning  to  the  receiver,  the  signal  is  said  to  double  bounce.  In  general, 
wave  pulses  that  deflect  off  of  an  odd  number  of  surfaces  are  called  odd  bounce,  and  wave 
pulses  that  deflect  off  an  even  number  of  surfaces  are  called  even  bounce  [10].  Physical 
boundary  conditions  on  perfect  conductors  lead  to  a  phase  shift  of  1 80°  on  every  bounce 
[9].  Therefore,  odd  and  even  bounce  scattering  are  determined  from  the  phase  difference 
of  the  co-polarization  terms,  S  hh  and  S  w  The  class  of  odd  bounce  scatterers  have  a 
scattering  matrix  of  the  form  [7,  9,  1 1] 


S 

—odd 


1 


0 

1 


(2.7) 


Backscatter  from  a  perfect  odd-bounce  conductor  leads  to  a  0°  shift  between  the 
polarization  channels,  since  the  co-polarization  terms  are  initially  offset  in  phase  by  180°. 
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Therefore,  the  co-polarization  scattering  magnitudes  are  equal,  |S  hh\  =  |5'  vvV  Similarly, 
the  class  of  even  bounce  scatterers  have  a  scattering  matrix  of  the  form  [7,  9,  1 1] 


S 

—even 


1 


1 

0 


0 

-1 


(2.8) 


Backscatter  from  a  perfect  even-bounce  conductor  leads  to  a  180°  phase  shift  between 
the  polarization  channels  [9],  since  two  180°  phase  shifts  cancel,  leaving  the  initial  180° 
shift.  Therefore,  the  co-polarization  scattering  magnitudes  are  opposite,  =  -|Syy|. 
Cross-polarization  scattering,  also  called  diffuse  scattering,  happens  when  the  reradiated 
pulse  is  orthogonal  to  the  one  transmitted  and  hence  bounce  parity  is  unknown.  The 
amount  of  diffuse  scattering  is  determined  by  the  cross-polarization  magnitude,  |S//y|.  The 
scattering  matrix  for  the  diffuse  class  of  scatterers  is  given  by  [7,  9,  11] 


S 

—diffuse 


0  1 
1  0 


(2.9) 


Common  canonical  shapes  have  a  known  scattering  response  and  are  often  used  for 
data  calibration  [10].  Canonical  scatterers  encompass  a  set  of  shapes  including  dihedrals, 
trihedrals,  spheres,  flat  plates,  cylinders  and  top-hats.  These  shapes  are  characterized  by 
several  parameters  including:  the  (v,  y,  z)  location  of  the  shape,  the  length  L,  the  height  H, 
the  radius  r,  and  the  roll,  pitch  and  yaw  angles  [10].  The  roll,  pitch,  and  yaw  angles 
define  the  counterclockwise  rotation  of  the  shape  with  respect  to  the  v,  y,  and  z  axis, 
respectively.  In  general,  dihedrals  and  top-hats  are  even  bounce  scatterers  and  plates, 
trihedrals,  cylinders,  and  spheres  are  odd  bounce  scatterers  [11].  Figure  2.1  shows  an 
example  of  all  these  shapes.  In  Chapter  4,  we  use  the  plate,  dihedral,  and  top-hat  to  illustrate 
the  proposed  polarization-preserving  image  enhancement  techniques. 

A  rectangular  plate  is  an  odd  bounce  scatterer  that  is  characterized  as  an  object 
with  one  flat  surface.  Reflected  energy  from  a  plate  is  directionally  dependent  [8,  10], 
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(b)  Dihedral  (c)  Trihedral 


Figure  2.1:  The  six  canonical  shapes  characterized  by  [10]. 


and  thus  backscattered  energy  is  only  received  by  a  monostatic  radar  when  the  plate  is 
normal  to  the  transmitted  signal.  A  plate  with  (roll, pitch, yaw)  =  (0°,0°,0°)  and  location 
(x,  y,  z)  =  (0, 0, 0)  is  as  pictured  in  Figure  2. la.  The  peak  RCS  cTpi^te  is  defined  as  [10] 

(Opiate  9  •>  (2.10) 


where  w  is  the  radar  wavelength. 

Dihedrals  are  even  bounce  scatterers  formed  from  two  flat  plates  joined  at  a  90°  angle. 
A  dihedral  with  (roll, pitch, yaw)  =  (0°,  0°,  0°)  and  location  (x,y,  z)  =  (0, 0, 0)  is  as  pictured 
in  Figure  2.1b.  The  peak  RCS  cr^ihedrai  is  given  by[10] 

(^dihedral  9  •  (2.11) 

The  top-hat  is  an  even  bounce  scatterer  that  is  formed  from  a  circular  flat  plate  and 
cylinder.  A  top-hat  with  (roll,pitch,yaw)  =  (0°,0°,0°)  and  location  {x,y,z)  =  (0,0,0)  is  as 
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pictured  in  Figure  2.1e.  The  top-hat  has  a  peak  RCS  (Ttophat  given  by  [10] 


^tophat  ~  /—  •  (2-12) 

w^f2 

2.2.2  Polarimetric  Decompositions  and  Classification. 

In  general,  a  scattering  matrix  S  corresponds  to  a  complex  coherent  target  that  is 
a  composite  of  many  scatterer  responses  [11].  Information  about  the  target  cannot  be 
obtained  purely  from  S.  Instead,  useful  information  is  obtained  by  decomposing  5  into  a 
combination  of  scattering  responses  of  canonical  shapes  [7,  11].  The  Pauli  decomposition 
is  one  standard  basis  used  to  decompose  a  scattering  matrix  [1,7].  It  characterizes  5  in 
terms  of  the  standard  scattering  matrices  of  a  sphere  5  ,  dihedral  S  ,  and  dihedral  with 

—odd  —even 

-45°  pitch  angle  5  .  Specifically,  the  Pauli  decomposition  is  given  as  [7,  9] 

— diffuse 

S  =  aS  +bS  +cS  (2.13) 

—  — odd  — even  — diffuse 

where 


a 

b 

c 


s  HH  Syy 


V2 

s  HH  -  Svv 

flS  HV- 


(2.14) 

(2.15) 

(2.16) 


The  Pauli  basis  decomposition  constants  are  such  that  |ap  is  the  power  scattered  by  targets 
with  odd  bounce,  is  the  power  scattered  by  targets  with  even  bounce  and  |cp  is  the 
power  from  diffuse  scattering  [7].  Similarly,  the  Krogager  decomposition  characterizes  the 
scattering  matrix  by  the  responses  of  a  sphere,  diplane  and  helix  [7]. 

The  Pauli  and  Krogager  polarimetric  decompositions  provide  a  simple  target 
classification  based  on  bounce  type.  Knowledge  of  the  bounce  type  is  useful  to  distinguish 
between  natural  and  man-made  targets  [12],  and  for  automatic  target  recognition  algorithms 
which  look  for  salient  features  to  classify  a  target  [11,  12].  Moreover,  knowledge  of  S 
allows  for  full  polarimetric  characterization  of  the  scattering  for  the  target  of  interest  [8], 
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(a)  (b)  (c) 

Figure  2.2:  The  (a)  gray-scale  SPAN  (b)  color  coded  SPAN,  and  (c)  Pauli  basis 
decomposition  of  a  dihedral. 

and  hence,  the  scattering  response  can  be  determined  in  any  polarimetric  basis.  Changing 
the  polarimetric  basis  allows  one  to  determine  the  scattered  polarization  state  for  any 
incident  polarization  state,  perhaps  to  enhance  specific  features  of  interest,  by  applying 
signal  processing  techniques  instead  of  physical  antenna  changes  [8,  9]. 

2.2.3  Polarimetric  Image  Formation. 

Since  a  fully  polarimetric  radar  is  a  multichannel  system,  all  the  polarimetric  channels 
must  be  considered  to  determine  the  total  scattered  power,  or  SPAN,  at  a  given  pixel.  The 
SPAN  is  given  by  [7] 

SPAN  =  IShh^  +  2\SHvf  +  l*^yyP.  (2.17) 

The  SPAN  can  be  presented  as  a  gray-scale  image,  where  the  intensity  of  a  given  pixel 
represents  the  relative  SPAN  magnitude  at  that  pixel.  Figure  2.2(a)  shows  an  example 
of  this  type  of  image.  While  this  method  does  effectively  image  a  scene  of  interest,  the 
polarimetric  information  is  not  portrayed. 

The  SPAN  can  also  be  presented  as  a  color-coded  image  where  the  relative  term 
intensities  represent  the  color  weights  on  a  Red,  Green,  Blue  (RGB)  scale  [7,  8].  For 
example,  let  15 hh?'  correspond  to  red,  2|5 hv?'  correspond  to  green,  and  |S  vv\^  correspond 
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to  blue.  Figure  2.2(b)  shows  an  example  of  this  type  of  image.  The  polarimetric  information 
is  portrayed  in  this  image;  however,  the  physical  interpretation  of  the  colors  can  be  difficult 
to  comprehend.  Furthermore,  notice  that  in  both  gray-scale  and  color-coded  SPAN  images, 
no  information  is  portrayed  about  the  relative  phase. 

One  method  to  present  both  the  polarimetric  information  and  relative  phase  is  through 
a  scattering  matrix  decomposition.  For  instance,  the  Pauli  decomposition  can  be  presented 
as  a  color-coded  image  where  |ap,  \b^,  and  |cp  represent  the  relative  weights  on  an  RGB 
scale  [1,  7].  In  this  case,  the  polarimetric  information  is  portrayed  with  a  meaningful 
physical  interpretation  and  the  relative  phase  is  captured  within  a  and  b.  Notice  that  SPAN 
=  lap  +  \b^  +  |cp,  so  the  total  scattered  power  is  maintained.  Figure  2.2(c)  shows  an  example 
of  a  dihedral  in  the  Pauli  basis  decomposition.  Notice  that  the  dihedral  appears  green,  since 
it  is  an  even  bounce  structure. 


2.3  Synthetic  Aperture  Radar 

A  synthetic  aperture  radar  (SAR)  is  a  radar  system  that  simulates  a  large  antenna 
aperture  by  moving  a  small  antenna  over  a  large  aperture  extent  [4].  The  dual  nature  of  SAR 
allows  for  creation  of  high  resolution  ground  images.  Images  are  created  of  the  measured 
ground  reflectivity  of  a  scene  of  interest.  There  are  two  SAR  modes,  stripmap  and  spotlight, 
which  distinguish  the  scene  geometry  based  on  the  radar.  In  stripmap  mode,  the  radar  beam 
points  in  a  fixed  direction  for  the  desired  aperture  duration.  In  spotlight  mode,  the  radar 
beam  points  at  a  fixed  target  for  the  desired  aperture  duration  [4]  [2].  Circular  flight  paths 
are  optimal  for  spotlight  mode  since  data  is  collected  for  a  full  360°  aperture.  This  work 
assumes  the  use  of  spotlight  mode  SAR  with  a  known  radar  flight  path. 

Typically,  a  SAR  transmits  a  linear  frequency  modulated  (LFM)  chirp  s{t)  given  by 


[4,  13] 


s{t) 


0, 


else 


(2.18) 
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Figure  2.3:  Illustration  of  grazing  angle  9  and  azimuth  angle  0. 

where  /o  is  the  carrier  frequency,  a  is  the  chirp  rate,  and  is  the  time  duration  of  the  chirp. 
The  received  signal  r(t)  is  first  mixed  with  delayed  in-phase  (I)  and  quadrature  (Q)  versions 
of  the  transmitted  chirp  and  then  low-pass  filtered.  Denote  this  signal  rdt).  Thus,  [4,  13] 

rdt)  -  W  f  (2.19) 

^  J-Ui 

where  A  is  a  scale  factor  that  accounts  for  attenuation  effects,  g(u)  is  the  scene  reflectivity 
function,  u  is  the  distance  along  the  range  direction,  c  is  the  speed  of  light,  tq  is  the  time 
delay  from  the  radar  platform  to  the  scene  center,  and  the  interval  denotes  the 

radar  range  based  on  chirp  duration.  Notice  that  by  ignoring  the  quadratic  phase  skew  term 
the  desired  scene  reflectivity  function  g(u)  is  the  inverse  Fourier  transform  of  rc(0[13]. 
The  signal  rdt)  is  also  known  as  phase  history  data. 

Phase  history  data  is  measured  at  each  grazing  angle  6  and  azimuth  angle  0,  as  shown 
in  Figure  2.3,  where  (0,0,0)  corresponds  to  the  scene  center  [4].  Assume  the  grazing 
angle  is  constant  throughout  the  flight  path.  Backprojection  [14]  or  Polar  Reformat  [4] 
algorithms  are  used  to  accurately  combine  the  information  from  all  azimuth  samples  to 
create  an  image  of  the  scene.  The  backprojection  algorithm  is  based  on  tomographic 
imaging  principles  [14].  For  each  (x,y)  coordinate  in  the  scene  extent,  the  inverse  Radon 
transform  is  computed  to  sum  the  reflectivity  contributions  from  each  azimuth  angle  at 
that  point  [14].  Alternatively,  the  Polar  Reformat  algorithm  (PFA)  interpolates  the  range 
and  azimuth  samples  into  a  Cartesian  raster  before  performing  a  two-dimensional  Fourier 
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transform  [4].  This  research  uses  PFA  to  find  scene  reflectivity.  The  range  resolution  pr  and 
cross-range  resolution  pc  are  determined  by  the  bandwidth  and  aperture  extent.  Specifically, 


[4] 


Pr 


C 

2Br 


w 

2A0’ 


(2.20) 


where  Be  is  the  chirp  bandwidth,  w  is  the  radar  wavelength,  and  A<p  is  the  aperture  extent. 

The  ideal  phase  history  data  can  also  be  thought  of  as  a  complex  vector  r  of  signal 
magnitude  and  phase.  Then,  the  sampled  scene  reflectivity  x  is  the  Fourier  Transform  of  r. 
Hence, 


r  =  ^. 


(2.21) 


where  T  denotes  the  Fourier  Transform  matrix.  In  other  words,  r  is  a  linear  combination 
of  the  basis  T  with  coefficients  x.  The  complex  measured  phase  history  y  is  then  given  by 
[15] 

y  =  ^  =  O'Vx  =  (2.22) 

where  O  is  the  measurement  matrix  which  accounts  for  the  limited  nature  of  data  collection 
that  transforms  r  to  y  and  0  =  OT.  The  matrix  represents  the  radar  point  spread 

function,  which  indicates  the  multi-dimensional  impulse  response  of  the  imaging  system. 
The  observed  scene  reflectivity  T“'y  is  a  convolution  of  the  point  spread  function  with  x. 
Ideally,  the  scene  reflectivity  x  is  recovered  by  deconvolution  of  the  point  spread  function. 


2.4  The  Assumption  of  Sparsity 

In  reality,  the  received  signal  contains  a  large  amount  of  interference  due  to  clutter, 
thermal  noise,  and  smearing  from  the  impulse  response.  Due  to  bandlimited  and  aperture- 
limited  data  collection,  0  is  ill-posed  and  therefore  not  invertible  [3].  It  follows  that 
directly  using  scattering  matrix  decomposition  theory  on  the  received  signal  produces 
images  that  may  give  incorrect  information  to  the  user,  such  as  a  misclassified  target.  The 
smearing  elfects  first  need  to  be  reduced  from  the  received  signal.  If  it  is  assumed  that 
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the  scene  of  interest  has  discrete  scatterers,  then  the  received  signal  will  be  sparse  and  can 
be  represented  in  a  basis  with  only  a  few  nonzero  entries  [3,  15].  Sparse  imaging  is  a 
regularization  method  that  estimates  the  scene  reflectivity  by  putting  signal  energy  in  only 
high  amplitude  pixels,  thus  emphasizing  targets  and  reducing  noise  and  smearing. 

The  assumption  of  sparsity  holds  when  the  measurement  matrix  O  satisfies  the 
Restricted  Isometry  Property  (RIP)  [3].  RIP  describes  matrices  which  are  orthogonal  or 
close  to  orthogonal.  Specifically,  O  should  satisfy  the  relationship  [3] 

(i-4)iwi  <  \m\i  <  (i+4)iwi,  (2.23) 

where  6^  is  the  restricted  isometry  constant  and  x  is  the  vector  of  measured  data.  The  radar 
collection  model  has  sufficient  performance  guarantees  which  validate  the  use  of  sparse 
reconstruction  [3]. 

The  system  of  equations  in  Eq.  (2.22)  is  under-determined.  The  sparse  estimate  for 
the  measured  data,  x,  is  given  by 

X  =  argmin  ||x||o  s.t.  ||y  -  ^||2  <  ei,  (2.24) 

X  —  - 

where  the  £2  norm  preserves  the  accuracy  of  the  signal  and  €i  is  a  predetermined 
infinitesimal  scalar  corresponding  to  the  level  of  error  allowed  [3]  [15].  However,  the  Co 
norm  is  an  NP-complete  problem  and  is  therefore  impractical  for  real-world  applications. 
Switching  the  4  to  an  Ci  norm  finds  a  relaxed  sparse  solution  [15].  By  including 
the  constraint  directly  in  the  objective  function,  the  sparsity  can  be  enforced  through  a 
regularization  parameter,  A.  This  produces  the  least  squares  problem  [3,  15] 

X  =  argmin  llj  -  +  dlWIi.  (2.25) 

X  —  - 

Note  that  the  Lagrange  multiplier  T  is  a  positive  scalar  [16].  The  ideal  value  of  A  for  a 
particular  scenario  can  be  found  though  a  dual  optimization  problem  [17].  This  research 
considers  sparse  enhancement  over  various  values  of  A. 
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2.5  Multichannel  Enhancement 


Often  times,  data  for  a  scene  of  interest  is  collected  over  multiple  channels,  such 
as  frequency,  elevation,  or  polarization.  Since  the  channels  carry  common  and/or 
complementary  information  about  the  same  scene,  combining  the  channels  proves  useful 
for  improving  estimates  of  appropriate  scene  attributes  [3].  For  example,  data  from 
multiple  elevations  can  be  used  to  refine  height  estimates  of  scatterers  within  a  scene 
[5].  Multichannel  images  are  formed  using  signal  processing  techniques  to  enhance  and 
combine  multiple  image  channels.  Enhancement  is  typically  done  independently  for 
each  channel.  Recent  research  has  also  considered  joint  enhancement,  which  involves 
simultaneously  enhancing  each  channel  based  on  a  known  interchannel  relationship  [5,  6]. 
In  [5],  scatterer  height  estimates  are  shown  to  be  more  accurate  using  joint  enhancement 
versus  independent  enhancement  techniques.  Furthermore,  in  [6]  the  interchannel 
covariance  matrix  is  used  to  improve  signal  estimates  and  reduce  noise  in  hyperspectral 
imaging.  Our  work  specifically  considers  image  enhancement  through  the  assumption  of 
sparsity. 

2. 5. 1  Independent  Enhancement. 

Independent  enhancement  is  appropriate  for  a  single  image  or  multiple  independent 
images.  In  the  case  of  polarimetric  data,  independently  enhancing  each  polarimetric 
channel  may  corrupt  the  original  signal  or  phase  ratio  between  the  channels  by  modifying 
energy  placement  in  the  pixels.  The  outcome  is  an  inferior  solution  that  ignores  the 
interchannel  coupling. 

Sparse  solution  estimates  for  n  multichannel  images  each  with  N  pixels  are  found 
independently  for  each  channel  through  the  least  squares  problem  [18] 

n 

=  arg  min  V  C(x.),  (2.26) 

where  the  objective  function  C(x.)  is  given  by 

C(x.)  =  ||y.-0x,]i  +  T,]|x.||i.  (2.27) 
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Hence,  a  particular  enhanced  image  x  based  on  the  phase  history  y  is  given  by 


x  =  argmin  ||y  -  0x11^  +  /l||jr||i. 

X  —  - 


(2.28) 


Most  numeric  solutions  of  Eq.  (2.21)  involve  determining  the  gradient  of  C.  To  avoid 
differentiating  the  £i  norm  for  values  of  x  close  to  zero,  ||x||i  is  approximated  as  [5,  18] 

N 

ll^lli  «  Yj^\x(j)\^  +  (2.29) 

7=1 


It  follows  that 


N 


X  ~  arg  min  (y  -  &x)^(y  -  &x)  +  A  ^(|^(7)P  +  e)*^^ 

■  ~  ~  ~  7=1 


N 


=  arg  min  y^y  -  2y^&x  +  ^^&x  +  A  ^(k(7)P  + 
■  ~  ~  ~  ~  7=1 

The  gradient  with  respect  to  x  is  given  by  [5] 


N 


Vx  =  -20^y  +  20^0jr  +  d  ^ 


l(^);l 


^  (|(X)/  +  6)1/2 


-20^j  +  20^0x  +  Ar(x)x, 


(2.30) 

(2.31) 


(2.32) 

(2.33) 


where  r(x)  =  Note  that  the  approximate  Hessian  matrix  is  given  by 

20^0  +  Tr(x)  [18]. 

The  least  squares  problem  in  Eq.  (2.26)  has  no  closed  form  solution  for  d  0 
[16].  However,  Eq.  (2.26)  can  also  be  phrased  as  a  linear  program  or  second  order 
cone  program  which  are  well  known  convex  optimization  problems  [16].  Optimization 
schemes  such  as  conjugate  gradient  descent  or  quasi-Newton  methods  can  be  used  to  find 
a  minimizing  solution  [16].  However,  these  methods  may  be  slow  to  converge  due  to 
expensive  computations.  In  [16],  an  algorithm  is  developed  using  Optimization  Transfer  to 
iteratively  converge  while  avoiding  expensive  computations.  The  output  of  each  iterate  k 
is  an  estimate  of  the  enhanced  image,  Within  the  (k  +  iterate,  a  surrogate  function 
M(x;  X*)  is  formed  and  minimized.  Hence, 


=  argminM(x;x*). 

X 


(2.34) 
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The  minimum  of  a  proper  surrogate  function  will  converge  to  the  minimum  of  the  original 
objective  function  [16].  In  other  words, 


2.5.2  Joint  Enhancement. 

Suppose  the  channels  are  known  to  be  coupled  through  a  set  of  functions.  In  this  case, 
it  is  appropriate  to  include  the  coupling  functions  as  additional  optimization  constraints. 
Current  joint  optimization  problems  are  formulated  as  [5]  [19] 

n 

X, ,  x, , . . . ,  X  =  arg  min  V  C (x  ) 

n 

s.t.  =  l...iV,  (2.38) 

i=\ 

where  the  real-valued  coupling  functions  are  continuous  and  N  is  the  number  of 

image  pixels.  Note  that  the  coupling  functions  have  an  independent  additive  form,  which 
implies  that  the  individual  images  may  be  still  be  independently  enhanced  [5]. 
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A  method  for  solution  involves  incorporating  the  constraint  functions  into  the 
objective  function  through  the  use  of  Lagrange  multipliers  [5].  Assume  that  a  local 
minimum  point  exists.  Hence,  there  exists  a  Lagrange  multiplier  yS'  for  each 

set  of  constraints  ^2/^2)’  •••’  which  corresponds  to  the  local  minimum.  It 

follows  that  the  analogous  unconstrained  objective  function  is  given  by  [5] 

n 

Y^C{x)+^h,  (2.39) 

where  h  is  the  stack  of  constraint  functions  given  by 

~  n  n  « 

h=  ,  (2.40) 

-  i=\  i=\  i=\ 

and is  a  vector  such  that  =  >0' .  Moreover,  (2.39)  is  locally  convex  around  (xj , ...,  ^). 

Using  the  Local  Duality  Theorem,  Eq.  (2.38)  is  equivalent  to  the  dual  function  [5] 

n 

x.,Xj,...,x  =  argmax  arg  min  >  C{x)  +  0^ h,  (2.41) 

for  any  yS  locally  near  yS'. 

Notice  that  the  optimization  problem  in  (2.41)  involves  solving  for  both  the  Lagrange 
multipliers  ySy  and  the  enhanced  images  x..  The  Dual  Descent  algorithm  [5]  is  a  numerical 
optimization  scheme  which  employs  a  dual  gradient  descent  to  solve  for  the  Lagrange 
multipliers  and  enhanced  images.  In  particular,  the  Dual  Descent  algorithm  fixes  the  set  of 
ySy  while  solving  for  the  restored  images  x.  and  then  subsequently  fixes  the  set  of  restored 
images  x.  to  solve  for  the  set  of  ySy  [5].  The  process  is  repeated  until  the  relative  change  in 

both  ySy  and  the  restored  images  has  been  sufficiently  minimized. 

Specifically,  the  estimated  Lagrange  multiplier  vector  yS  is  updated  for  each  iterate  k 
using  the  gradient  based  equation  [5] 

(2.42) 

where  is  the  set  of  constraint  functions  based  on  the  {k  -  1)'*  image  estimates, 
...,00^.  Subsequently,  the  enhanced  image  estimates  for  an  iterate  k  are  found  by 
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Input:  Measured  images  j  ,  forward  operators  0  ,  sparseness  parameters 


p,  Ai,  stepsize  a,  and  convergence  parameters  62,  where  i  =  1, n. 


Output:  Jointly  enhanced  images  x.,  where  /  =  1, n. 
Initialize:  =  y  ,0^  =  1,  and  =  Y!i=\ 


Algorithm: 
while  e  >  62  do 

Calculate  =0  +  aVi^k 
while  Relative  change  in  all  do 

Form  diagonal  matrix  F  using 
Solve  for  given  0*^ 

Calcuate  relative  change  in  given  as 
end  while 

Save  solution  as  =  ^0 
Calculate 

Calculate  relative  change  in  yS  given  as  e 


Ill'll 


Cost  fn  using  0^^-  Cost  fn  using  0 
Cost  fn  using  0 


end  while 


Output  x.  =  x^+* 


Figure  2.4:  Pseudocode  for  the  Dual  Descent  Algorithm  [5]. 


minimizing  (2.35)  using  conjugate  gradient  or  quasi-Newton  methods.  Table  2.1  shows  the 
pseudocode  for  the  Dual  Descent  method.  The  Dual  Descent  method  provides  a  starting 
point  to  solve  the  joint  enhancement  problem  for  polarimetric  radar  data. 

2.5.3  Issues. 

Recent  research  into  joint  enhancement  has  shown  promising  results  for  multichannel 
images  with  similar  support  regions,  such  as  signal  magnitude  or  phase  [5,6].  In  the  case  of 
polarimetric  data,  specific  scatterers  have  high  signal  strength  in  some  channels  but  not  in 
others.  This  is  true  for  a  dihedral,  which  has  similar  support  in  the  S hh  and  S  vv  channels. 
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(a)  (b) 

Figure  2.5:  The  (a)  15 hh\  (=  |5  vv\),  and  (b)  |5 hv\  unenhanced  responses  of  a  dihedral. 

but  has  a  low  response  in  the  S  hv  channel.  Figure  2.4  shows  the  scattering  responses  of  the 
dihedral. 

Moreover,  the  formulation  of  coupling  constraints  in  Equation  (2.38)  does  not  account 
for  interchannel  relationships  in  which  an  independent  additive  form  cannot  be  achieved. 
For  example,  consider  the  convex  function 

=  ]^\z^+ z^\.  (2.43) 

Notice  that  without  further  assumptions  (2.43)  cannot  be  written  in  the  form  h{z^,z^  = 
h\{z^  +  h2iZ2),  where  h\,h2  are  any  convex  functions.  Similarly,  the  relationship 
between  polarimetric  channels  given  by  scattering  matrix  decomposition  does  not  follow  an 
independent  additive  form.  Consequently,  the  generalized  joint  enhancement  theory  must 
be  modified  before  application  to  polarimetric  data.  Chapter  3  presents  such  a  generalized 
method. 
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111.  Joint  Enhancement  of  Polarimetric  Data 


RAdar  data  may  be  collected  over  multiple  polarization  channels  [3].  Channels  may 
carry  common  and/or  complementary  information  about  the  scene  of  interest. 
Combining  appropriate  channels  produces  improved  estimates  of  scene  attributes  [5]. 
Multichannel  images  are  formed  by  enhancing  and  combining  multiple  image  channels. 
This  research  specifically  considers  image  enhancement  using  the  assumption  of  sparsity. 

Although  image  enhancement  is  usually  done  independently  for  each  channel,  recent 
research  has  showed  promising  results  for  joint  image  enhancement  [5]  [6].  Joint  image 
enhancement  accounts  for  the  interchannel  relationship  that  exists  when  multiple  channels 
are  coupled  through  a  known  set  of  functions.  The  Dual  Descent  algorithm  is  a  current 
scheme  designed  to  optimize  the  tradeoff  between  sparsity  in  the  image  and  preservation 
of  interchannel  relationships  that  satisfy  Eq.  (2.38)  [5].  Joint  enhancement  has  proved 
successful  for  multichannel  datasets  with  similar  support  regions,  such  as  signal  magnitude 
or  phase  [5,  6,  19].  In  the  case  of  polarimetric  radar  data,  the  received  signal  may 
be  different  from  one  channel  to  the  next.  However,  the  polarimetric  channels  are 
related  to  one  another  and  together  provide  useful  information  about  the  scattering  of  a 
target  of  interest.  It  is  important  to  preserve  the  interchannel  relationship  during  image 
enhancement. 

This  section  derives  two  functional  relationships  to  describe  polarimetric  channel 
coupling  among  received  signal  channels.  Then,  an  optimization  problem  is  constructed 
to  maintain  signal  fidelity,  enforce  sparsity,  and  preserve  interchannel  coupling.  The 
convexity  of  the  solution  is  analyzed.  Finally,  an  algorithmic  solution  of  the  optimization 
problem  is  developed. 
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z 


Figure  3.1:  Illustration  of  a  polarization  state  vector  mapped  onto  the  Poincare  Sphere 
shown  in  red,  and  scaled  sub-sphere  shown  in  black. 

3.1  Non-Convex  Constraint  for  Polarimetric  Data 

Recall  that  a  fully  polarimetric  monostatic  radar  produces  three  distinct  polarimetric 
data  channels:  Shh,  Shv  and  Sw  For  clarity,  let  y^^  and  x^y,  Xyy 

represent  the  received  phase  history  and  enhanced  scene  reflectivity  of  the  S  hh,  S  hv  and 
Svv  channels,  respectively.  The  three  enhanced  polarimetric  channels  and 

x^,^,  are  the  foundation  for  a  variety  of  radar  applications  including  data  calibration  and 
target  classiflcation,  and  therefore,  need  to  be  accurate  signals  with  precise  interchannel 
relationships. 

The  Poincare  Polarization  sphere  [7,  9]  provides  a  useful  picture  of  all  possible 
scattered  polarization  states.  In  particular,  each  point  on  the  sphere  corresponds  to  a 
speciflc  scattered  polarization  state.  Recall  that  the  scattered  polarization  state  can  be 
fully  described  by  the  incident  polarization  state  and  the  scattering  matrix.  Therefore, 
preservation  of  the  scattering  matrix  is  important  to  maintain  the  proper  direction  vector  of 
the  point  on  the  Poincare  sphere,  i.e.,  the  precise  scattered  polarization  state.  Although 

the  radius  of  the  Poincare  sphere  is  informative,  the  direction  vector  fully  describes  the 
polarization  state  [7].  Permitting  the  radius  to  change  but  maintaining  the  vector  direction 
allows  for  signal  sparsity  without  changing  the  inherent  information  about  the  polarization 
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state.  To  preserve  the  vector  direction,  the  contributions  of  the  enhanced  channels  should 
remain  proportional  to  the  contributions  of  the  received  channels.  Moreover,  the  phase 
between  channels  is  informative  for  the  classification  of  scattering  types.  Since  phase 
is  captured  within  a  complex  exponential,  the  ratio  of  any  two  channels  gives  the  phase 
differences  as  well  as  the  ratio  of  the  signal  magnitudes.  Therefore,  we  propose  an 
appropriate  coupling  constraint 


—HV  —V\A 


^hhU) 

+ 

^HvU) 

+ 

XhvU) 

lyyii) 

XyyU) 

XyyU) 

XhhU) 

to  preserve  relative  channel  magnitudes  and  phase  differences.  The  real-valued  constraint 
of  Eq.  (3.1)  is  a  measure  of  the  absolute  error  between  received  and  enhanced  polarimetric 
ratios  for  each  resolution  cell  j,  or  pixel,  in  the  image.  Polarimetric  preservation  occurs 
when  x^y,  Xyy)  =  0  for  each  pixel  j. 

3.1.1  Convexity. 

Consider  the  convexity  of  hj.  Clearly,  hj  is  not  a  linear  function.  In  general,  a  function 
is  convex  if  the  Hessian  is  positive  semidefinite  [20].  The  Lagrangian 


is  given  by 


—HV’  —V\A 


dhj  dhj 


dhj  1 

dXyy{j)\ 


(3.2) 
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where 


\  ^w(j)  Zyy^-^^  ^V/y(7)  j 


(3.3) 

1  Zyy^-^^  y^yyU)XfJHU)  -  ^^J^j)Xyy{j)  X^y{j)  (^^^0)  ^ 

l^yy(7)l  l^yy(.^3l  (^yy(7)^yi/(7)  ~  ^^^(7)^yy(7))  _  ^hvO) 

(3.4) 


dhj  _  d  f 

_|_ 

Zffy^7)  x^y(;) 

_l_ 

Zffy^7)  x^y(;)  ] 

dXfjyU)  dxjjy(j)  1 

^yy(7)  ^yy(7) 

^yy(7)  Xyyij) 

(3.5) 


1  ^yyU)^HV^j^  ^^y(j)^yy(j)| 

l^yy(./)l  l^yy(.(^l  ^^yy^j^—HV^j^  ~  ^fiy^j^—VV^j^^ 

1  ^^^(i)^//y(./)  ~  ^fjyU)^HH^j^ 

I— ffff (.(^1  (Z////^33^//y(7)  ~  ^fjy(J)^HH^j^^ 

^  I  _l_  -HV^-^^  _  — gy(33  Zgy^-^^ 

^—VV^j^  1^  —VV^j^  ^YV^^^  —VV^j^ 

I _  -HH^j^  \  / _  &V^\ 

^  ^gg(;)  I  Ivv^j)  XyyU)  j  ^  Xjjy(j)  [lyyU)  Xyy(J)  j 

xIyU)  I  I  xIyU)  Izhv^ ^//vO)  I 

lyyU)  Xyy(j)  y^^(J)  Xyy(J) 


(3.6) 

(3.7) 

^/ry(7)  ] 
^hhU)  I 

(3.8) 

(3.9) 
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Therefore  the  Hessian 


^  is  given  by  ^  = 

— h  — h 


d^hj 

d^hj 


d^hj 


d^hj 


d^hi 


'5^hv0') 

d^hi 


where 


d% 

d^hi 


—  9v  (  A  (  _  -HV^j^  \ 

-HV^J^[y^^(j)  x^^U)) 
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dx^nU^dx^yU) 


■HH 

y 


U) 

O’) 


^HV 


U) 


( _  y-Hv^j^  \ 

-////O)  / 


\ 


d% 


dXfJH(j)dXyy(j) 

d% 


XhhU) 
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/Zzzzz' 

ij)  XfjnU)  ^ 

\^yy^j^  ^yy(7)  / 

J^^(7)  XfjuU) 

±vv^J' 

y_yy^j)  XyyU) 

dxly(j) 


=  0 


d^hj 


dXfjy{j)dXyy{j) 


-VV 


(j) 


uif) 


d% 

dxlyU) 


y.yyU) 


d^hi 


d^hj 


^—HV 

^yiyU) 
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vU) 
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Xyyij) 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


Let  /  =  [/"  f^  /f  ]"  €  Cl  Then  Hf  = 


a  b  c 
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-9v  (  A  (  _  -HV^i'>  \ 

-^vyj>\y_^^(j)  x„„(j)) 
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XfjyU) 

XhhU) 

^hhU)  \ 
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Xi) 


Xhh^J) 
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(3.17) 
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C  =  /l 


/ \ 


y*2  /  y\ 

±vvX]) 

y_vvd) 

y.yyU) 

+  /2 


vU) 


t/0) 


y*2  /  j\ 

^//y(7) 
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^yy(j) 

+  /s 


( 

O  Y* 
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(3.18) 
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(3.19) 
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Equation  (3.19)  simplifies  to 


fHf  =  /f/i 


'H  ^ 


-9v  (  A  i  _  -ffyO)  \ 


9^y0*)  ^//vO) 
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1 

At  this  point,  it  becomes  difficult  to  make  generalizations  about  the  convexity  of  the 
function  based  on  the  definiteness  of  the  Hessian.  However,  notice  that  hj  does  have 
discontinuities  whenever  any  components  in  x.  are  zero,  i  e  {HH,  VV}.  Therefore,  the 
domain  of  each  image  would  certainly  need  restriction  in  some  fashion  to  be  convex. 
However,  an  argument  can  be  made  for  local  convexity  of  the  function  during  optimization 
if  an  initial  guess  for  each  x.  is  sufficiently  close  to  the  desired  relative  minimum  point  x*. 
Such  an  intital  guess  may  be  the  received  image  y ,. 


3.2  Convex  Constraint  for  Polarimetric  Data 

Equivalently,  polarization  preservation  occurs  when  the  received  phase  history 
polarization  ratios  are  equal  to  the  enhanced  scene  reflectivity  polarization  ratios,  i.e.. 
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U) 


y^^U) 


^HvU) 


y_HvU) 


^yy(2)  ^yy(2) 


IhvU) 


(3.21) 
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Cross-multiplication  leads  to  another  appropriate  real- valued  constraint  given  by 


8i 


-h 


y |^^yO)^yyO) 


(3.22) 


Notice  that  each  term  within  g —ny^  -vy^  ^  composition  of  two  convex  functions: 
the  absolute  value  function  and  a  linear  function  of  the  channels.  Convexity  is  preserved 
under  composition  and  addition,  so  it  follows  that  gj  is  convex  [20].  To  see  this  rigorously, 
consider  the  Hessian  of  gjiXj^ff,  x^^y,  Xyy).  The  Lagrangian  ^gjiXfffj,  x^^y,  Xyy)  is 

^ gji^HH-'  —HV^  —VV^  ~ 


ijy'H.HH'’  —HV'  "yy^ 

dgj  dgj  dgj 


dx^v^j)  dXyy(j)\ 


where 


^gj^^HH'  —HV’  —VV^ 

dxHnU) 


(3.23) 


(3.24) 


^gj(^HH^^HV’^VV^  ,  Zyy^-^^-UvO) 


dx^vU) 


(3.25) 


^gj(^HH^^HV’^VV^  _  y  y,yyU)XHHU) 

dXyyU)  y^^U)Xyy(j)-y^Jj)XfjHU) 


-VV 


(3.26) 


w^Jiyyvj/  :LyyW/:2i//yv 

Notice  that  the  fractional  term  of  each  partial  derivative  simplifies  to  +1,  meaning  each 
partial  derivative  is  not  explicitly  a  function  of  the  corresponding  partial  term.  Therefore 
the  Hessian  H  is  given  by 

— g 

0  0  0 

0  0  0.  (3.27) 

0  0  0 


H  = 

— g 
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Clearly,  f^Hf  =  0  for  all  /  e  which  implies  the  Hessian  is  positive  semidefinite. 


3.3  Polarimetric  Data 

The  joint  image  enhancement  problem  for  polarimetric  radar  data  can  be  phrased  as 


mm 
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y  -  0x„„ 

-HH  ^^HH 
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-HV  ^^HV 
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y  -  0x,,,, 

—VV  - k  k 

^  Q'2 

subject  to  (3.28) 

WXyyWl  <  QTS 

P A—hh'  —hv  —  T 

where  for  small  values  of  or  1,0^2, 0^3,  sparse  solutions  are  enforced  and  for  small  values  of 
T,  interchannel  relationships  are  preserved  through  the  coupling  constraint  pj  =  {hj,  gj}  for 
each  pixel  j.  Note  that  the  last  constraint  in  Eq.  (3.28)  could  be  an  equality  constraint, 
equal  to  zero,  which  would  preserve  interchannel  coupling.  However,  the  requirement  is 
relaxed  to  allow  for  small  error  tolerance  since  the  signals  will  be  noisy.  Moreover,  notice 
that  the  same  0  and  A  are  assumed  for  each  channel  since  all  the  channels  are  collected  at 
the  same  time  and  therefore  undergo  the  same  data  collection  effects. 

Using  the  method  of  Lagrange  multipliers  and  assuming  the  Kuhn-Tucker  conditions 
hold  [20],  the  optimal  solution  satisfies 


VL  =  0, 


(3.29) 


where 


L  -  L\  +^ipi(x*  ,  X*  ,  X*)  +  fi2P2(x*f,fr,  X*,  X*  )  H - I3nPn{^„„,  x*,  Xyy),  (3.30) 
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(3.31) 


x*^H^  x^y,  Xyy  are  the  relative  minimum  points,  A  is  the  Lagrange  multiplier  corresponding 
to  the  ^i-norm  constraint,  and  ySj  is  the  Lagrange  multiplier  corresponding  to  the 
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polarization-preserving  constraint  Xyy)  for  the  jth  pixel.  Since  there  are  =  3 

data  channels  and  N  pixels  per  channel,  the  gradient  of  the  Lagrangian  becomes 
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(3.32) 


PnVLhh’’  —hv  —v'\A 
where  we  have  approximated  the  (’i-norm  as  in  [18]. 

Using  the  Local  Duality  Theorem,  Eq.  (3.28)  can  be  converted  to  the  equivalent  dual 
problem 


^vv  =  min  Li (x^^,  x^y,  Xyy)  +  0^p,  (3.33) 


where  fils  a  vector  such  that  fi{j)  =  Pj  and  pQ)  =  pj. 
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3.4  Algorithmic  Solution 

To  solve  the  polarimetric  joint  enhancement  problem  given  in  Eq.  (3.33),  we  employ 
an  iterative  dual  gradient  method  similar  to  the  Dual  Descent  method  described  in  Table 
(2.1)  [5].  Moreover,  we  set  the  initial  estimates  of  x.  as  the  corresponding  y  .  Other  inputs 
include  the  desired  sparsity  weighting  A,  data  collection  operator  0,  stepsize  a,  and 
other  desired  convergence  parameters.  Then,  a  double  loop  is  initiated,  where  we  alternate 
solving  for  the  estimates  x.  and  the  Lagrange  multiplier  estimates  fi. 

The  image  estimates  x.  remain  fixed  while  solving  for  fi.  Specifically,  the  fi  update 
equation  for  the  Ath  iterate  is 

^  (3.34) 

where  is  the  set  of  constraint  functions  based  on  the  (k  -  ly^  image  estimates, 

^Hv'^vv-  Then  each  entry  yS/  in  is  scaled  by  the  received  signal  strength  sj  given 
by 

Sj  =  +  yjyv^j))^  (3.35) 

in  order  to  emphasize  preservation  on  regions  of  the  scene  where  received  signal  strength  is 

^  Jc  1 

larger.  Next,  the  remain  constant  while  solving  for  the  image  estimates  Specifically, 
a  preconditioned  conjugate  gradient  method  with  a  quadratic  surrogate  function  is  used, 
as  in  [16].  Chapter  4  demonstrates  the  suitability  of  this  algorithm,  considering  both  the 
non-convex  constraint  pj  =  hj  and  the  convex  constraint  pj  =  gj. 
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IV.  Results 


This  research  considers  optimization  constraints  to  preserve  coupling  between 
multiple  polarization  channels  of  radar  data  collected  for  the  same  scene  of  interest.  For 
a  given  resolution  cell  or  pixel  j,  scattering  types  are  preserved  through  the  non-convex 
constraint 

h(r  r  r  .  _  X^yU) 

^  ^ 

or  the  convex  constraint 

8ji^HH^—HV^—VV^  ~  VvOl  ~  +  J^yO)^yyO)  “ 

~  •  (4-2) 

In  this  section,  the  joint  optimization  problem  in  Section  3.4  is  solved  using  the  algorithm 
described  in  Section  3.5  for  both  the  convex  and  non-convex  coupling  constraints. 
Polarimetric  SAR  datasets  of  several  canonical  scatterers  [10]  are  considered.  Jointly 
enhanced  images  are  compared  with  corresponding  images  found  through  independent 
enhancement.  Specifically,  the  objective  function  is  broken  down  and  individual  costs  are 
considered  to  demonstrate  fidelity  of  signal  estimates  as  well  as  preservation  of  polarimetric 
channel  coupling.  Moreover,  a  scattering  measure  of  success  is  derived  and  compared  for 
the  enhanced  images. 

4.1  Measures  of  Error 

Recall  that  the  optimization  objective  function  L  is  given  by 

2  2  2 

^  =  Ihh~  2  +  Ihv~  2  +  Zyy  -  ^vv  ^  +  MXj,y\h  +  As\\Xyy\h 

l^\P\ ^^hh’  —hv’  ^yy)  l^2P2i8LHH'’  —HV’  — yy^  -I-  •  •  •  -I-  ^yy)?  (4.3) 

where  Pjix^fj,  x^y,  Xyy)  is  either  the  convex  constraint  gj  or  the  non-convex  constraint  hj. 
The  individual  costs  of  sparsity,  error  of  signal  estimates,  and  error  of  polarimetric  preser- 
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vation  are  considered  and  are  determined  from  the  respective  terms  of  L.  Additionally,  the 
error  of  scattering  classification  is  derived  from  the  Pauli  basis  decomposition. 

4.1.1  Objective  Function. 

Signal  fidelity  indicates  the  accuracy  of  the  signal  estimate.  The  error  due  to  loss  of 
signal  fidelity  for  a  particular  estimate  x.,  i  e  {HH,  HV,  VV}  is  given  by  the  2-norm  cost 

(4.4) 

The  amount  of  sparsity  present  in  a  signal  estimate  x.  is  given  by  the  1-norm  cost 

11^, 111.  (4.5) 

The  error  due  to  loss  of  polarimetric  preservation  is  calculated  using  the  cost  of  the 
polarimetric  coupling  constraint.  Hence,  when  hj  was  enforced  during  optimization,  the 
error  cost  is  given  by  x^^y,  Xyy),  and  similarly  for  gj. 

4.1.2  Scattering  Measure. 

Preservation  of  the  scattering  matrix  is  important,  particularly  for  target  classification 
efforts.  Along  with  demonstrating  that  the  objective  function  is  effectively  enforced,  the 
classification  of  the  enhanced  images  is  considered.  Forming  an  image  from  the  Pauli  basis 
decomposition  creates  a  RBG  vector  c  =  [R  G  B]^  for  each  pixel,  where  high  values  of 
red  correspond  to  strong  odd  bounce  reflection,  high  values  of  green  correspond  to  strong 
even  bounce  reflection,  and  strong  values  of  blue  correspond  to  high  orthogonal  reflection. 
If  little  to  no  returned  signal  is  detected,  c  has  values  close  to  0  and  then  the  pixel  is  black. 
For  clarity,  off-target  black  pixels  are  recolored  white.  The  classification  error  at  a  given 
pixel  is  given  by 

\\c-c\\l,  (4.6) 

where  c  is  the  RGB  vector  of  the  enhanced  signal  and  c  is  the  RBG  vector  of  the  ideal 
scattering  matrix  at  that  pixel. 
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4.1.3  Parameters. 


This  research  considers  the  measures  of  error  over  several  parameters,  including 
sparsity  level,  noise  level,  and  specific  regions  of  the  scene.  In  each  case,  error  is  averaged 
over  500  Monte  Carlo  noise  simulations. 

The  value  of  the  sparsity  enforcing  parameter  A  ranges  from  0.1  -  1.3.  As  d  is 
increased,  the  level  of  sparsity  enforced  is  increased.  In  [17],  an  algorithm  is  presented 
for  determining  the  value  of  A  corresponding  to  the  lowest  objective  function  cost.  For  this 
research,  an  appropriate  range  of  A  values  was  picked  visually. 

Furthermore,  both  datasets  are  simulated  and  therefore  do  not  contain  noise.  White 
Gaussian  noise  (WGN)  is  added  to  the  data  at  peak  signal-to-noise  Ratio  (SNR)  ranging 
from  20  dB  to  50  dB,  with  20  dB  indicating  a  high  noise  environment  and  50  dB  indicating 
a  low  noise  environment. 

Moreover,  measures  of  error  were  considered  over  specific  regions  of  the  scene.  The 
full  scene  error  is  the  error  over  the  entire  image.  The  on-target  error  is  the  error  over 
image  masks  defined  where  the  normalized  span  is  at  least  0.5  dB.  Similarly,  the  off-target 
error  is  the  error  over  image  masks  defined  where  the  normalized  span  is  lower  than  0.5 
dB,  corresponding  to  all  pixels  not  on-target.  The  maximum  pixel  error  is  the  error  at  the 
pixel  with  the  strongest  scattered  power. 

4.2  Canonical  Shapes 

A  complex  target  can  be  broken  down  into  the  known  scattering  responses  of  simple 
shapes,  or  canonical  scatterers,  as  discussed  in  Section  2.2.  To  illustrate  polarimetric 
preservation,  we  consider  scenes  containing  odd,  even,  and  diffuse  bounce  canonical 
shapes.  Shape  data  is  simulated  over  an  azimuth  extent  of  20°,  elevation  of  30°,  wavelength 
of  0.03m,  center  frequency  of  fc  =  10e9  Hz  and  bandwidth  B  =  2.95e9  Hz. 
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Figure  4.1:  Top  Row:  The  configuration  for  the  plate,  dihedral,  titled  dihedral,  and  group 
scenes.  Second  Row:  The  SPAN  of  the  plate,  dihedral,  tilted  dihedral  and  group  scenes, 
along  with  regions  of  interest. 


Polarization  preservation  imaging  results  are  found  for  an  individual  plate,  an 
individual  standard  dihedral,  an  individual  tilted  dihedral,  and  a  group  scene  consisting 
of  a  plate,  standard  dihedral  and  top-hat.  For  clarity,  the  standard  dihedral  scene  is  called 
the  ‘dihedral  scene’,  and  the  tilted  dihedral  scene  is  called  the  ‘tilted  dihedral  scene.’  To 
eliminate  comparisons  across  RCS,  the  shape  parameters  are  chosen  such  that  the  RCS  is 
the  same  for  each  canonical  shape.  Figure  4.1  shows  the  scene  setup  and  boxes  in  the 
on-target  mask  (black)  and  maximum  pixel  (star)  for  each  scene.  The  off-target  mask 


Table  4.1:  Canonical  Shape  Scene  Parameters. 


Parameters 

Plate 

Dihedral 

Top-hat 

Tilt-Dihedral 

(x,y,z) 

(0, 0, 0)  m 

(0, 0, 0)  m 

(-l,-l,-1.5)m 

(0, 0, 0)  m 

L 

0.8m 

0.6m 

- 

0.6m 

H 

0.4243m 

0.4m 

1.6478m 

0.4m 

r 

- 

- 

Im 

- 

RCS 

1608.5  m2 

1608.5  m2 

1608.5  m2 

1608.5  m2 

Roll/Pitch/Yaw 

30°  pitch 

- 

- 

45°  pitch 
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contains  all  pixels  outside  the  on-target  mask.  Table  4.1  lists  the  shape  parameters  for 
each  canonical  scene.  It  should  be  noted  that  the  tilted  dihedral  has  the  same  parameters  as 
the  dihedral  along  with  a  45°  pitch  angle.  Moreover,  the  group  scene  consists  of  a  top-hat 
with  parameters  as  listed  as  well  as  a  plate  at  (1, 1, 2)  m  and  standard  dihedral  at  (0,0, 2) 
m.  All  other  plate  and  dihedral  parameters  are  the  same  for  the  group  scene. 


Figures  4.2  -  4.5  show  the  |5 hh\,  |5' hv\  and  |S  vv\  signal  responses  for  each  scene  at 
each  stage  of  processing  using  30  dB  of  WGN  and  X  =  0.7.  The  pure  simulated  images 
in  the  top  row  of  each  figure  clearly  show  the  impulse  response  smearing.  Notice  that  due 
to  the  orientation  of  all  the  scatterers,  tilted  dihedral  is  the  only  object  with  significant 
response  in  the  HV  channel.  The  second  row  of  each  figure  shows  the  images  after 
adding  noise  at  a  moderate  SNR  of  30  dB.  The  third  row  in  Figures  4.2-4.5  shows  the 
results  of  independently  enhancing  each  polarization  channel  under  a  sparsity  constraint 
with  X  =  0.7.  The  fourth  row  shows  the  results  of  jointly  enhancing  the  polarimetric 
images  using  the  non-convex  constraint  hj.  Finally,  the  fifth  row  shows  the  results  of  joint 
enhancement  using  the  convex  constraint  gj. 

First  consider  noise  level.  Notice  that  in  general  the  jointly  restored  images  (both 
convex  and  non-convex)  have  reduced  the  noise  better  than  the  independently  restored 
images;  in  the  titled  dihedral  scene  each  enhancement  type  removes  noise  well.  Next, 
consider  signal  restoration.  In  general,  the  independent  images  and  non-convex  jointly 
enhanced  images  produce  similar  results  for  the  HH  and  VV  channels;  the  convex  jointly 
enhanced  images  appear  to  have  preserved  some  high-level  smearing  close  to  the  target, 
while  decreasing  the  strength  of  the  signal  in  the  on-target  region;  this  phenomena  is  most 
pronounced  in  the  single  shape  cases,  and  less  dramatic  in  the  group  scene.  Moreover,  in 
the  dihedral  case  the  non-convex  constraint  appears  to  have  successfully  restored  some  of 
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Figure  4.2:  Plate  scene  images  for  all  stages  of  enhancement.  Top  Row:  Simulated  \S hhI 
\S Hvl  and  |5  vv\  with  no  added  noise.  Second  Row:  Simulated  |5 hh\,  I-S hvI  and  \S  vv\ 
with  30  dB  of  WGN.  Third  Row:  Independently  restored  IS hhI  |S hv\,  and  |S  vv\  with  no 
polarimetric  coupling  constraint  for  A  =  0.7.  Fourth  Row:  Jointly  restored  |5 hh\,  |S hv\, 
and  |5  vv\  with  non-convex  coupling  constraint  for  A  =  0.7.  Last  Row:  Jointly  restored 
\S hhI  \S hv\,  and  \S  vv\  with  convex  coupling  constraint  for  A  =  0.7. 
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Figure  4.3:  Dihedral  scene  images  for  all  stages  of  enhancement. Top  Row:  Simulated 
\S hhI  I-S Hvl  and  |5  vv\  with  no  added  noise.  Second  Row:  Simulated  \S hhI  I-S hvI  and 
|5  vv\  with  30  dB  of  WGN.  Third  Row:  Independently  restored  IS  hhI  |S  hvI  and  |5  vv\ 
with  no  polarimetric  coupling  constraint  for  A  =  0.7.  Fourth  Row:  Jointly  restored  |5 hh\, 
\S Hvl  and  |5  vv\  with  non-convex  coupling  constraint  for  A  =  0.7.  Last  Row:  Jointly 
restored  \S hhI  \S hv\,  and  \S  vv\  with  convex  coupling  constraint  for  A  =  0.7. 
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Figure  4.4:  Tilted  Dihedral  scene  images  for  all  stages  of  enhancement.Top  Row: 
Simulated  and  \Svv\  with  no  added  noise.  Second  Row:  Simulated 

and  |5yy|  with  30  dB  of  WGN.  Third  Row:  Independently  restored  |5f/y|, 

and  |5  yy|  with  no  polarimetric  coupling  constraint  for  A  =  0.7.  Fourth  Row:  Jointly 
restored  |5 hh\,  I-S hvI  and  \S  yy|  with  non-convex  coupling  constraint  for  A  =  0.7.  Last 
Row:  Jointly  restored  |5 hh\,  I-S hvI  and  |5  yy|  with  convex  coupling  constraint  for  A  =  0.7. 
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Figure  4.5:  Group  scene  images  for  all  stages  of  enhancement.Top  Row:  Simulated  \S hhI 
\S Hvl  and  |5  vv\  with  no  added  noise.  Second  Row:  Simulated  |5 hh\,  I-S' hvI  and  \S  vv\ 
with  30  dB  of  WGN.  Third  Row:  Independently  restored  IS hhI  |S hv\,  and  |S  vv\  with  no 
polarimetric  coupling  constraint  for  A  =  0.7.  Fourth  Row:  Jointly  restored  |5 hh\,  |S hv\, 
and  |5  vv\  with  non-convex  coupling  constraint  for  A  =  0.7.  Last  Row:  Jointly  restored 
\S hhI  \S hv\,  and  \S  vv\  with  convex  coupling  constraint  for  A  =  0.7. 
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Figure  4.6:  The  2-norm  cost  for  the  plate  scene  as  a  function  of  A  (SNR  =  30  dB).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  A  (SNR  =  30  dB). 


the  HV  channel.  However,  for  the  group  case  the  non-convex  constraint  appears  to  have 
added  some  HV  response  not  present  in  the  original  images,  which  is  undesirable.  The 
tilted  dihedral  has  good  HH,  HV,  and  VV  response  for  all  three  enhancement  types. 

4.2.1  Signal  Fidelity. 

The  error  due  to  loss  of  signal  fidelity  is  considered  for  each  canonical  scene  over 
all  the  target  regions  for  increasing  values  of  A  and  SNR.  Figure  4.6  compares  the  2- 
norm  costs  of  the  independently  restored  plate  images  with  both  the  non-convex  jointly 
restored  plate  images  and  the  convex  jointly  restored  plate  images  for  increasing  values 
of  A.  Similarly,  Figure  4.7  compares  the  2-norm  costs  of  the  independently  restored  plate 
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Figure  4.7:  The  2-norm  cost  for  the  plate  scene  as  a  function  of  SNR  (A  =  0.7).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  SNR  (A  =  0.7). 


images  with  the  non-convex  jointly  restored  plate  images  and  the  convex  jointly  restored 
plate  images  for  increasing  SNR.  Notice  that  the  2-norm  cost  is  considered  individually  for 
the  HH,  HV  and  VV  channels.  Moreover,  note  that  in  the  graph  legend  Ind  represents  an 
independently  enhanced  channel,  Joint  C  represents  a  convex  jointly  restored  channel,  and 
Joint  NC  represents  a  non-convex  jointly  restored  channel.  The  subgraphs  do  not  have  the 
same  scale.  Figures  4.8-4.9,  4.10-4.11,  and  4.12-4.13  show  the  corresponding  images  for 
the  dihedral,  tilted  dihedral  and  group  scenes,  respectively. 

First  consider  the  independently  restored  images.  For  each  scene,  the  independently 
restored  images  have  a  consistently  lower  full  scene  error  cost  in  each  channel  compared  to 
the  jointly  restored  image  channels.  Furthermore,  the  on-target  error  of  the  independently 
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Figure  4.8:  The  2-norm  cost  for  the  dihedral  scene  as  a  function  of  A  (SNR  =  30  dB).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Off-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  A  (SNR  =  30  dB). 


enhanced  channels  is  consistently  lower  than  the  corresponding  off-target  error.  In  the 
single  shape  scenes,  the  three  independently  restored  channels  have  a  similar  on-target  error 
magnitude.  The  three  independently  restored  channels  display  on-target  error  magnitude 
differences  when  the  number  of  scatterers  increases.  Specifically,  the  group  scene  displays 
a  lower  independently  restored  HV  channel  error  compared  to  the  HH  and  VV  channels, 
which  is  attributable  to  the  one  region  of  low  scattering  response  in  the  HV  channel 
compared  to  the  distinct  regions  of  high  scattering  response  in  the  HH  and  VV  channels. 
Conversely,  the  tilted  dihedral  scene  displays  a  higher  full  scene  error  in  the  independently 
restored  HV  channel  versus  the  HH  or  VV.  Upon  inspection,  the  higher  error  in  the  HV 
channel  comes  from  the  off-target  region  and  not  the  on-target  region.  In  general,  the 
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Figure  4.9:  The  2-norm  cost  for  the  dihedral  scene  as  a  function  of  SNR  (A  =  0.7).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  SNR  (A  =  0.7). 


2-norm  cost  steadily  increases  in  each  channel  as  A  increases  and  steadily  decreases  in 
each  channel  as  SNR  increases.  These  trends  make  sense;  the  signal  fidelity  and  sparsity 
have  a  tradeoff.  Thus  higher  penalties  for  sparsity  lead  to  an  increased  error  in  signal 
fidelity.  Moreover,  as  additive  noise  decreases  the  signal  becomes  easier  to  restore  leading 
to  decreased  error  in  signal  fidelity. 

Next  consider  the  non-convex  jointly  restored  images.  The  non-convex  jointly  restored 
image  channels  consistently  show  the  HV  channel  with  the  most  full  scene  error,  followed 
by  the  HH  channel  and  then  VV  channel.  For  the  plate  and  dihedral  scenes,  the  three  non- 
convex  jointly  restored  channels  have  approximately  the  same  error;  in  the  tilted  dihedral 
and  group  scene  there  is  some  difference  between  the  HV  channel  and  the  other  two 
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Figure  4. 10:  The  2-norm  cost  for  the  tilted  dihedral  scene  as  a  function  of  A  (SNR  =  30  dB). 
The  (Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target 
and  (Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  A  (SNR  =  30  dB). 

channels.  Moreover,  notice  that  the  on-target  error  is  significantly  less  than  the  olf-target 

error,  which  is  desirable.  In  all  scenes,  the  on-target  non-convex  jointly  restored  channels 
have  a  similar  level  magnitude  of  error  to  the  independent  channels,  particularly  for  small 
A.  Moreover,  notice  that  for  the  plate,  dihedral  and  group  scenes,  the  non-convex  jointly 
restored  channels  have  less  on-target  error  than  the  corresponding  convex  jointly  restored 
channels.  However,  in  the  tilted  dihedral  scenes  the  non-convex  jointly  restored  channels 
have  a  distinctly  higher  on-target  error  than  the  corresponding  convex  channels.  Next, 
notice  that  the  non-convex  jointly  restored  channels  converge  to  the  independently  restored 
channels  as  SNR  increases.  This  is  desirable,  since  the  independent  channels  minimize 
error  in  a  no-noise  environment,  so  it  follows  that  the  non-convex  jointly  restored  channels 
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Figure  4.11:  The  2-norm  cost  for  the  tilted  dihedral  scene  as  a  function  of  SNR  (A  =  0.7). 
The  (Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target 
and  (Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  SNR  (T  =  0.7). 


should  converge  to  a  similar  solution.  Overall,  the  non-convex  jointly  restored  images 
display  similar  signal  accuracy  compared  to  the  independent  case. 

Finally,  consider  the  convex  jointly  restored  image  channels.  The  convex  jointly 
restored  channels  consistently  show  the  HH  and  VV  channels  with  the  same  error  and 
the  HV  channel  with  consistently  lower  error.  Moreover,  notice  that  in  the  plate  and 
dihedral  scenes,  the  full  scene  error  of  the  convex  jointly  restored  HH  and  VV  channels 
is  significantly  higher  than  the  corresponding  full  scene  error  in  the  independent  and  non- 
convex  jointly  restored  channels.  Looking  closer,  a  significant  portion  of  the  convex  jointly 
restored  full  scene  error  is  due  to  on-target  error,  which  is  significantly  higher  than  its 
off-target  error.  This  trend  is  undesirable;  the  on-target  regions  of  the  scene  are  arguably 


46 


Figure  4.12:  The  2-norm  cost  for  the  group  scene  as  a  function  of  A  (SNR  =  30  dB).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Off-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  A  (SNR  =  30  dB). 

more  important  than  the  off-target  regions.  In  the  tilted  dihedral  scene,  the  convex  jointly 
restored  channels  perform  better  in  the  on-target  region  than  the  off-target.  Finally,  notice 
that  while  all  error  types  of  the  convex  jointly  restored  channels  do  decrease  with  SNR,  we 
do  not  see  convergence  to  the  corresponding  independent  channels. 

Overall,  the  independently  restored  channels  have  a  consistently  lower  2-norm  cost 
than  the  corresponding  non-convex  jointly  restored  channels  and  convex  jointly  restored 
channels.  In  the  scenes  with  strong  odd  and  even  bounce  scatterers,  the  HV  channel  has 
consistently  lower  2-norm  cost  than  the  NH  and  VV  channels  for  each  enhancement  type. 
Moreover,  the  non-convex  jointly  restored  channels  have  a  lower  2-norm  cost  compared 
to  the  convex  jointly  restored  channels  in  the  scenes  with  strong  odd  and  even  bounce 
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Figure  4.13:  The  2-norm  cost  for  the  group  scene  as  a  function  of  SNR  (A  =  0.7).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Off-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  2-Norm  for  increasing  SNR  (A  =  0.7). 


scattering.  However,  in  the  scenes  with  strong  diffuse  scattering  the  ffV  channel  tends 
to  have  more  error  than  the  corresponding  ffff  and  VV  channels,  and  we  see  a  lower  2- 
norm  cost  in  the  convex  jointly  restored  channels  compared  to  the  independently  restored 
channels. 

42.2  Sparsity. 

The  cost  of  the  sparsity  term  is  considered  for  each  scene  over  increasing  values 
of  A  and  SNR  for  various  scene  regions.  Figure  4.14  compares  the  1-norm  cost  of  the 
independently  restored  plate  images  with  both  the  non-convex  jointly  restored  plate  images 
and  convex  jointly  restored  plate  images  for  increasing  values  of  A  and  various  scene 
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Figure  4.14:  The  1-norm  cost  for  the  plate  scene  as  a  function  of  A  (SNR  =  30  dB).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  A  (SNR  =  30  dB). 


regions.  Similarly,  Figure  4.15  compares  the  1-norm  cost  of  the  independently  restored 
plate  images  with  both  the  non-convex  jointly  restored  plate  images  and  convex  jointly 
restored  plate  images  for  increasing  SNR  over  various  scene  regions.  Notice  that  the  1- 
norm  cost  is  considered  individually  for  the  NH,  HV,  and  VV  channels.  Figures  4.16-4.17, 
4.18-4.19,  and  4.20-4.21  show  the  corresponding  images  for  the  dihedral,  tilted  dihedral 
and  group  scenes,  respectively. 

First  consider  the  independently  restored  channels.  For  each  scene,  and  each  scene 
region,  the  independently  restored  images  have  a  higher  1-norm  cost  in  each  channel 
compared  to  the  corresponding  joint  convex  and  joint  non-convex  channels.  The  exception 
is  the  olf-target  cost  in  the  plate  and  dihedral  scenes,  where  the  convex  jointly  restored 
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Figure  4.15:  The  1-norm  cost  for  the  plate  scene  as  a  function  of  SNR  (A  =  0.7).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  SNR  (A  =  0.7). 


images  have  a  higher  cost  than  the  independent  images.  From  Figures  4.2  -  4.5,  it  follows 
that  the  higher  costs  are  due  in  part  to  the  higher  noise  level  in  the  independent  images 
compared  to  the  corresponding  joint  images,  and  in  part  to  a  stronger  on-target  signal. 
Furthermore,  notice  in  plate,  dihedral  and  group  scenes,  the  independently  restored  images 
have  a  significantly  higher  signal  response  in  the  ffff  and  VV  channels  compared  to  the  ffV 
channel,  which  leads  to  a  consistently  lower  1-norm  cost  for  the  NV  channel  compared  to 
the  Nff  and  VV  channels  in  all  scene  regions.  In  the  tilted  dihedral  scene,  where  the 
NV  channel  has  a  strong  signal  response,  the  independently  restored  NV  channel  has  a 
higher  1-norm  cost  than  the  NN  and  VV  channels.  Considering  the  average  error  per  pixel 
indicates  the  independently  restored  images  have  a  low  cost  average  per  pixel  in  the  off- 
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Figure  4.16:  The  1-norm  cost  for  the  dihedral  scene  as  a  function  of  A  (SNR  =  30  dB).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Off-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  A  (SNR  =  30  dB). 


target  region,  and  a  high  average  cost  per  pixel  in  the  on-target  region.  Hence,  the  sparse 
pixels  are  picked  primarily  from  the  on-target  region,  which  is  desirable.  Finally,  notice 
that  the  1-norm  cost  of  the  independent  channels  steadily  decreases  over  A  and  decreases 
over  SNR,  as  expected. 

Next  consider  the  non-convex  jointly  restored  images.  The  non-convex  jointly  restored 
channels  display  similar  trends  to  the  independent  case;  specifically,  in  the  plate,  dihedral 
and  group  scenes,  the  HH  and  VV  channels  have  a  consistently  higher  sparsity  cost  than 
the  HV  channel;  in  the  tilted  dihedral  scene,  the  HV  channel  has  a  consistently  higher 
cost.  Moreover,  the  non-convex  jointly  restored  channels  converge  to  the  sparsity  level  of 
the  corresponding  independent  channels  as  SNR  increases.  Based  on  Figures  4.2-4.5,  the 
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Figure  4.17:  The  1-norm  cost  for  the  dihedral  scene  as  a  function  of  SNR  (A  =  0.7).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Off-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  SNR  (A  =  0.7). 


independent  and  convex  jointly  restored  images  converge  to  a  similar  solution.  Next,  notice 
that  the  non-convex  jointly  restored  images  have  a  strong  on- target  signal  and  a  low  signal 
elsewhere,  which  manifests  as  a  high  average  pixel  1-norm  cost  in  the  on-target  region  and 
a  low  average  pixel  1-norm  cost  in  the  off  target  region,  as  desired. 

Finally,  consider  the  convex  jointly  restored  images.  The  convex  jointly  restored 
channels  display  similar  trends  to  the  independent  case;  specifically,  in  the  plate,  dihedral 
and  group  scenes,  the  HH  and  VV  channels  have  a  consistently  higher  sparsity  cost  than 
the  HV  channel;  in  the  tilted  dihedral  scene,  the  HV  channel  has  a  consistently  higher 
cost.  In  the  plate,  dihedral,  and  group  scenes,  the  convex  jointly  restored  channels  have  a 
lower  1-norm  cost  average  per  pixel  in  the  on-target  region  compared  to  the  corresponding 
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Figure  4.18:  The  1-norm  cost  for  the  tilted  dihedral  scene  as  a  function  of  A  (SNR  =  30  dB). 
The  (Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target 
and  (Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  A  (SNR  =  30  dB). 


independent  and  non-convex  jointly  restored  channels,  and  a  higher  1-norm  cost  average 
per  pixel  in  the  off-target  region.  From  Figures  4.2-4.5,  its  clear  that  the  convex  jointly 
restored  images  have  a  lower  signal  strength  on-target,  and  have  preserved  some  impulse 
response  smearing  off-target.  In  the  titled  dihedral  scene,  the  convex  jointly  restored 
channels  have  a  higher  pixel  average  1-norm  cost  on-target  compared  to  the  non-convex 
jointly  restored  channels. 
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Figure  4.19:  The  1-norm  cost  for  the  tilted  dihedral  scene  as  a  function  of  SNR  (A  =  0.7). 
The  (Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Olf-Target 
and  (Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  SNR  (T  =  0.7). 


54 


Figure  4.20:  The  1-norm  cost  for  the  group  scene  as  a  function  of  A  (SNR  =  30  dB).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Off-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  A  (SNR  =  30  dB). 
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Figure  4.21:  The  1-norm  cost  for  the  group  scene  as  a  function  of  SNR  (A  =  0.7).  The 
(Top  Left:)  Full  scene  error  (Top  Right:)  On-Target  Error  (Lower  Left:)  Off-Target  and 
(Lower  Right:)  Max  Pixel  cost  of  the  1-Norm  for  increasing  SNR  (A  =  0.7). 


56 


4.2.3  Polarimetric  Preservation. 


The  error  due  to  loss  of  polarimetric  preservation  is  considered  for  each  canonical 
shape  over  increasing  values  of  A  and  SNR  for  several  scene  regions.  Figure  4.22  compares 
the  polarimetric  constraint  cost  of  the  independently  restored  plate  images  with  both  the 
non-convex  jointly  restored  plate  images  and  convex  jointly  restored  plate  images  for 
increasing  values  of  A  and  various  scene  regions.  Similarly,  Figure  4.23  compares  the 
polarimetric  constraint  cost  of  the  independently  restored  plate  images  with  both  the  non- 
convex  jointly  restored  plate  images  and  convex  jointly  restored  plate  images  for  increasing 


SNR  over  various  scene  regions.  Notice  that  the  polarimetric  constraint  cost  is  considered 


individually  for  the  terms 
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4.25,  4.26-4.27,  and  4.28-4.29  show  the  corresponding  images  for  the  dihedral,  tilted 


dihedral  and  group  scenes,  respectively.  Notice  that  among  all  figures,  similar  scattering 


in  the  HH  and  VV  responses  manifests  itself  as  a  similar  error  for  the  ratios 
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First  consider  the  independently  restored  image  ratios.  The  independently  restored 
polarimetric  ratios  have  a  high  error  compared  to  the  corresponding  jointly  restored  ratios 
for  all  scenes  and  scene  regions.  In  the  tilted  dihedral  scene,  all  preservation  types  appear 
to  have  similar  levels  of  error  in  all  scene  regions.  This  improvement  in  the  independently 
restored  channels  may  be  due  to  the  similar  signal  strengths  in  each  channel.  In  all  scenes, 
the  independently  restored  ratios  have  a  high  full  scene  error,  the  bulk  of  which  comes  from 
the  olf-target  region.  However,  the  independently  restored  on-target  polarimetric  ratio  error 
is  high  compared  to  the  corresponding  jointly  restored  on- target  polarimetric  ratio  error. 
For  the  plate  and  dihedral  scenes,  the  independently  restored  image  ratio  error  tends  to 
increase  as  A  increases  and  decreases  as  SNR  increases,  as  expected.  However,  for  the 
tilted  dihedral  and  group  scenes,  the  full  scene  and  off-target  error  of  the  independently 
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Figure  4.22:  The  polarimetric  preservation  cost  for  the  plate  scene  as  a  function  of  A  (SNR 
=  30  dB).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Olf-Target 
Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent  and 
jointly  enhanced  images  for  increasing  values  of  A. 


restored  image  ratios  decreases  as  A  increases.  Moreover,  notice  that  in  the  plate  and 
dihedral  scenes,  the  full  scene  and  off-target  error  of  the  independently  restored  ratios  and 
convex  jointly  restored  ratios  has  a  peak  at  30  dB,  which  may  be  due  in  part  to  the  large 
amount  of  30  dB  impulse  response  smearing  throughout  the  plate  and  dihedral  response. 

Next  consider  the  jointly  restored  image  ratios.  For  all  scenes,  the  non-convex  jointly 
restored  polarimetric  ratios  have  a  consistently  low  off  and  on-target  error  magnitude.  For 
all  scenes,  the  maximum  pixel  error  of  the  non-convex  jointly  restored  ratios  is  lower  than 
the  corresponding  independently  restored  ratios,  but  higher  than  the  corresponding  convex 
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Figure  4.23:  The  polarimetric  preservation  cost  for  the  plate  scene  as  a  function  of  SNR 
(A  =  0.7).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Olf-Target 
Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent  and 
jointly  enhanced  images  for  increasing  values  of  SNR. 


jointly  restored  ratios.  These  trends  are  desirable;  polarimetric  preservation  does  better  in 
the  jointly  restored  channels  compared  to  the  independently  restored  channels. 
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Figure  4.24:  The  polarimetric  preservation  cost  for  the  dihedral  scene  as  a  function  of  A 
(SNR  =  30  dB).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Off- 
Target  Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent 
and  jointly  enhanced  images  for  increasing  values  of  A. 
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Figure  4.25:  The  polarimetric  preservation  cost  for  the  dihedral  scene  as  a  function  of  SNR 
(A  =  0.7).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Olf-Target 
Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent  and 
jointly  enhanced  images  for  increasing  values  of  SNR. 
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Figure  4.26:  The  polarimetric  preservation  cost  for  the  tilted  dihedral  scene  as  a  function  of 
A  (SNR  =  30  dB).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Olf- 
Target  Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent 
and  jointly  enhanced  images  for  increasing  values  of  A. 
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Figure  4.27:  The  polarimetric  preservation  cost  for  the  tilted  dihedral  scene  as  a  function  of 
SNR  (i  =  0.7).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Olf- 
Target  Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent 
and  jointly  enhanced  images  for  increasing  values  of  SNR. 
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Figure  4.28:  The  polarimetric  preservation  cost  for  the  group  scene  as  a  function  of  A  (SNR 
=  30  dB).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Olf-Target 
Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent  and 
jointly  enhanced  images  for  increasing  values  of  A. 
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Figure  4.29:  The  polarimetric  preservation  cost  for  the  group  scene  as  a  function  of  SNR 
(A  =  0.7).  Top  Left:  Full  scene  error  Top  Right:  On-Target  error  Lower  Left:  Olf-Target 
Error  Lower  Right:  Max  Pixel  Error  of  polarization  preservation  of  the  independent  and 
jointly  enhanced  images  for  increasing  values  of  SNR. 
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4.2.4  Scattering  Measure. 

The  scattering  measure  is  considered  for  each  scene  region  over  increasing  values  of 
A  and  SNR  as  computed  by  Eq.  (4.6).  The  error  due  to  improper  scattering  classification 
is  considered  for  each  canonical  shape  over  increasing  values  of  A  and  SNR  for  the  on- 
target  region.  Figure  4.30  shows  the  Pauli  basis  decomposition  of  the  ideal  plate  on-target 
scattering  matrix,  independently  restored  images,  non-convex  jointly  restored  images  and 
the  convex  jointly  restored  images  for  A  =  0.7  and  SNR  =  30  dB.  Figure  4.31  compares 
the  plate  scattering  classification  error  of  the  independently  restored  plate  images  with  both 
the  non-convex  jointly  restored  plate  images  and  convex  jointly  restored  plate  images  for 
increasing  A  and  SNR.  Figures  4.32-4.33, 4.34-4.35,  and  4.36-4.37  show  the  corresponding 
images  for  the  dihedral,  tilted  dihedral  and  group  scenes,  respectively. 

First  consider  the  independently  restored  images.  The  independently  restored  images 
have  a  consistently  lower  scattering  classification  error  than  the  non-convex  jointly  restored 
images.  However,  notice  the  graph  scale  indicates  that  the  error  cost  difference  between 
the  all  the  restored  images  is  less  than  0.1.  The  scattering  error  of  the  independently 
restored  images  increases  as  A  increases,  and  decreases  as  SNR  increases  as  expected. 
The  independent,  convex  jointly  restored,  and  non-convex  jointly  restored  scattering 
classification  error  tend  to  converge  as  SNR  increases.  Notice,  however,  that  the  jointly 
restored  images  do  perform  better  in  low  SNR  cases  compared  to  the  independently 
restored  images.  In  particular,  the  convex  jointly  restored  images  have  a  uniform  scattering 
classification  error  across  A  and  SNR. 

4.3  Summary 

In  theory,  jointly  optimizing  polarimetric  radar  data  to  enforce  sparsity  and  preserve 
scattering  matrix  relationships  among  image  channels  should  lead  to  improved  target 
classification.  The  figures  and  graphs  throughout  Ch.  4  demonstrate  that  the  improved 
preservation  in  the  jointly  enhanced  channels  compared  to  the  independently  enhanced 


66 


channels  did  not  consistently  lead  to  improved  scattering  classification  elforts  under 
the  metric  considered.  This  discrepancy  may  arise  from  definition  differences  between 
the  polarimetric  coupling  constraint,  which  preserved  scattering  matrix  relationships  of 
the  received  signals,  and  the  scattering  classification  metric,  which  compared  the  ideal 
classification  to  the  estimated  classification.  However,  the  jointly  enhanced  images  did 
perform  better  in  scattering  classification  than  the  independently  restored  images  for  some 
low  SNR  cases.  These  cases  should  be  considered  further,  as  polarimetric  preservation 
may  provide  helpful  insight  into  a  scene  of  interest  in  high  noise  cases.  Moreover,  future 
research  could  consider  ways  to  rephrase  the  polarimetric  coupling  constraint  to  better 
reflect  imperfections  in  the  received  signals. 
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Figure  4.30:  Pauli  basis  decomposition  for  the  plate  scene.  The  Top  Left:  ideal  scattering 
matrix,  Top  Right:  independently  enhanced  images,  Lower  Left:  non-convex  jointly 
enhanced  images,  and  Lower  Right:  convex  jointly  enhanced  images  at  SNR  =  30  dB 
and  A  =  0.7. 


Figure  4.31:  On-target  average  scattering  error  for  the  plate  scene  over  A  and  SNR. 
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Figure  4.32:  Pauli  basis  decomposition  for  the  dihedral  scene.  The  Top  Left:  ideal 
scattering  matrix,  Top  Right:  independently  enhanced  images,  Lower  Left:  non-convex 
jointly  enhanced  images,  and  Lower  Right:  convex  jointly  enhanced  images  at  SNR  =  30 
dB  and  A  =  0.7. 


Figure  4.33:  On-target  average  scattering  error  for  the  dihedral  scene  over  A  and  SNR. 
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Figure  4.34:  Pauli  basis  decomposition  for  the  tilted  dihedral  scene.  The  Top  Left:  ideal 
scattering  matrix,  Top  Right:  independently  enhanced  images,  Lower  Left:  non-convex 
jointly  enhanced  images,  and  Lower  Right:  convex  jointly  enhanced  images  at  SNR  =  30 
dB  and  A  =  0.7. 


Figure  4.35:  On-target  average  scattering  error  for  the  tilted  dihedral  over  A  and  SNR. 
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Figure  4.36:  Pauli  basis  decomposition  for  the  group  scene.  The  Top  Left:  ideal  scattering 
matrix,  Top  Right:  independently  enhanced  images,  Lower  Left:  non-convex  jointly 
enhanced  images,  and  Lower  Right:  convex  jointly  enhanced  images  at  SNR  =  30  dB 
and  A  =  0.7. 


Figure  4.37:  On-target  average  scattering  error  for  the  group  scene  over  A  and  SNR. 
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V.  Conclusion 


Synthetic  aperture  radar  images  contain  impulse  response  effects  that  appear  to  smear 
or  blur  the  image  [4],  Sparse  imaging  is  a  regularization  method  that  provides  a  way 
to  overcome  limited  datasets  for  image  reconstruction  [3].  Sparsity  constraints  enforced 
alongside  data  fidelity  constraints  lead  to  signal  energy  placement  in  only  high  amplitude 
pixels  [3],  thus  emphasizing  high  magnitude  scatterers  in  the  scene  while  reducing  the 
effects  of  noise  and  the  impulse  response. 

Image  data  for  a  scene  of  interest  may  be  collected  over  multiple  polarization  channels 
[5].  In  the  case  of  polarimetric  SAR  images,  regularization  techniques  are  typically  applied 
independently  to  each  polarimetric  channel.  However,  independent  processing  does  not 
account  for  cross-channel  coupling  and  may  corrupt  the  polarimetric  information  in  the 
signals.  Recent  research  into  joint  enhancement  techniques  has  shown  promising  results 
for  multi-channel  datasets  with  similar  regions  of  signal  magnitude  and/or  phase  [5,  6]. 
However,  in  the  case  of  polarimetric  SAR  data,  scattering  may  be  present  in  some  channels 
and  not  in  others.  This  thesis  mathematically  formulated  multi-channel  sparse  imaging 
for  polarimetric  radar  data  using  a  joint  enhancement  algorithm  to  enforce  sparsity  and 
polarimetric  coupling  constraints. 

5.1  Summary  of  Findings 

Two  candidate  functional  relationships  were  separately  considered  to  describe 
polarimetric  coupling  among  received  signal  channels:  one  convex  function  g  and  one 
non-convex  function  h.  Specifically,  these  functions  were  designed  to  preserve  the  Sinclair 
scattering  matrix.  The  two  functions  g  and  h  were  reformed  as  optimization  constraints 
gj  and  hj,  respectively.  Then,  the  coupling  constraints  were  considered  separately  in 
an  optimization  problem  constructed  to  maintain  signal  fidelity,  enforced  sparsity,  and 
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preserve  inter-channel  coupling.  The  constrained  problem  was  turned  into  an  unconstrained 
dual  problem  similar  to  [5]  by  adding  the  sparsity  constraint  and  polarimetric  coupling 
constraint  to  the  objective  function  using  Lagrange  multipliers. 

An  iterative  dual  gradient  descent  algorithm  [5,  16]  was  used  to  alternatively 
calculate  the  updated  scene  estimates  for  each  channel  and  then  the  maximizing  Lagrange 
multipliers  for  each  coupling  constraint.  Results  were  found  for  several  polarimetric  SAR 
datasets,  including  simulated  radar  data  of  various  canonical  scatterers.  Jointly  enhanced 
images  were  compared  to  corresponding  images  found  through  independent  enhancement, 
taking  into  consideration  signal  fidelity,  sparsity,  polarimetric  preservation,  and  scattering 
classification. 

Overall,  the  jointly  enhanced  image  channels  displayed  significantly  better  polarimet¬ 
ric  preservation  compared  to  the  corresponding  independently  restored  image  channels. 
More  research  is  needed  to  understand  how  polarimetric  preservation  can  be  used  to  im¬ 
prove  target  classification. 

5.2  Future  Research 

Future  research  into  this  topic  could  consider  the  scattering  classification  improvement 
of  jointly  enhanced  polarimetric  images  using  more  extensive  classification  tests.  For 
example,  one  could  test  whether  dictionary  estimates  of  scatterers  are  improved  with 
jointly  enhanced  polarimetric  images  versus  independently  enhanced  polarimetric  images. 
Moreover,  future  research  could  consider  further  polarimetric  coupling  constraints 
designed  to  emphasize  desired  polarization  responses  in  the  image,  such  as  odd  bouce 
or  dihedral  reflection.  This  should  not  be  confused  with  enhancement  of  image  regions, 
edges,  or  points,  which  has  been  well  documented. 
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